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ABSTRACT 

We calculate the local UV flux density in the GALEX MIS FUV and NUV bands 
using redshifts provided by SDSS DR7. Luminosity functions are calculated for the 
overlapping MIS and SDSS sample, allowing flux densities to be measured and the 
local star formation rate (SFR) to be calculated using volumes much larger than 
previous FUV based estimates. We calculate flux densities for a dust corrected low 
redshift (0.013 ^ z ^ 0.1) sample of fu-FUV = 22.24 ± 3.13 xlO^^h ergs s"^ Hz"^ 
Mpc-^, fu-NUV = 38.54 ± 5.30 xlO'^^h ergs s"^ Rz'^ Mpc"^. The star formation 
rate density found is 0.0312 ± 0.0045 h M0yr~^ Mpc~^. This is larger than published 
rates recently found using the UV implied SFR, though the major discrepancy is the 
correction made for dust attenuation and once this is dealt with consistently the results 
agree well. These values are also consistent with recent Hq, derived SFRs. Once cosmic 
variance is taken into account most of the recent SFRs at low redshift {z ^ 0.3) found 
in the literature can be brought into agreement, however the lowest redshift values 
{z ^ 0.045) do appear to be significantly lower. 

Key words: GALEX: SDSS: Luminosity Function: Star Formation Rate 



1 INTRODUCTION 

In its simplest form the UV to near-IR output of the en- 
tire Universe can be predicted with remarkably few pre- 
scriptions. Primarily we need to know the initial mass 
function (IMF) for stars, be this constan t (e.g. IWvse 



1997 , and references therein) or varying fe.g. IWilkins et al 



20081): the evolution of s t ars once they ha ve formed (e.g 



Bruzual CharlotI l2003l : iMarastonI l2005l ) and the rate 



at which baryons have been turned into stars through- 
out the history of the Universe to our current epoch 
([Hopkins Beacoml 2006). On top of this we must un- 
derstand the role of dust in attenuating the light pumped 
into the cosmos (e .g. ICalzetti et aIll2Q0Ql : iTuffs et al.ll2QQ4l : 
[Driver et al.l l2008h . Though we can state this quite suc- 
cinctly, each component of the puzzle is challenging for both 
theory and observation. In our ideal picture we would be 
able to state that the Universe produces stars at a rate de- 
termined from the star formation history (SFH) and with 
a mass distribution given by the IMF. To predict the lo- 
cal Cosmic Energy Distribution (CSED), all we need to do 
is model the evolution of stars using a population synthe- 
sis model and account for the role of dust in absorbing and 
re-emitting the light in the mid and far-IR. 
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In practice, our best current observations and theoreti- 
cal models cannot perfectly measure or predict any of these 
separate components, and thus the z — ^ light output of the 
Universe is non-trivial to explain. In terms of predictions, 
weak constraints on IMFs hav e been ma de through high 
resolution modelling (jBonnell et al.l l200lh . Equally, popu- 
lation synthesis models have become increasingly sophisti- 
cated, but no purely theoretical models have created a star 
as we observe them, and major issues exist with the lat- 
ter stages of stellar evolution where dredge up and mass 
loss have a significant effe ct on the radiation output of stars 
but are difficult to model (|Marigo et al. 1 120081 ). The star for- 
mation rate ( SFR) has been predicted in a variety of sim- 
ulations (e.g. Cole et al 



been p re 
1 l20QQl : [f 



ilernquist k Springel 2003'; 

[Na gamine et al. 2 0061 ) but these are ultimately limited by 
our observations (because models are designed to fit the 
data) and the assumed IMF — clearly the conversion be- 
tween UV or Hq; output and SFR can differ between reality 
and simulation. These theoretical difficulties highlight the 
full complexity of the problem, and ultimately mean that 
observations are still the leading light in building our un- 
derstanding of the flux output of the Universe. 

Observationally the IMF can be measured to some ex- 
tent in star forming regions. This is difficult to do directly 
since the most massive stars are short-lived and have a de- 
layed proto- stellar period compared to less massive stars. 
This short window of existence makes it unlikely that they 
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will be observed directly, placing significant uncertainty on 
the high mass slope and turnover of the IMF. Star forming 
regions are also expected to be highly obscured by dust and 
much star forming activity will be hidden entirely. Despite 
these restrictions the IMF has been rneasur ed directly in var- 
ious regions of the Milky Way (|Scaldl 19981 ) , with the general 
finding thait it is not radically dissimilar to a Salpeter IMF 
("Salpeter 1955). Observations also play a key role in pro- 
ducing stellar models since spectral measurements of stars 
are used as templates in synthesis codes to evolve stars and 
galaxies. A vast amount of effort has been invested into mea- 
suring the cosmic SFH, surveys have been c onducted at a 
large range of redshifts (see co mpendiums bv lWilkins et al.] 
I2OO8I : iHopkins BeacQmll2006 . and references therein) util- 
ising different tracers in order to produce a SFR density 
as a function of look-back t ime, c reating the, now ubiqui- 
tous, Lilly- Mad au diagram (jLillv et al. 1996; Madau ~et al.l 
Il996l . ll998V Despite these tremendous efforts, the SFH of 
the Universe is not well constrained. In particular, the sig- 
nificant effects of cosmic variance (see Driver & Robotham 
2010) and even determining the true volume of space sur- 
veyed at high redshift makes constraining the SFR at early 
epochs extremely difficult. Current results allow the cosmic 
SFR to be roughly constant for z > 2 or even declining 
sharply, the single consensus is that it was not any higher in 
the earliest epochs of the Universe. As well as being a key 
ingredient in calculating the CSED, some form of IMF is re- 
quired to convert all star formation tracers into conventional 
units of SFR. Often this conversion is done implicitly, but 
there is enough convergence in different IMF models that 
the resultant SFRs are not hugely sensitive to the IMF even 
for different tracers. 

The vast majority of UV output in the FUV and NUV 
bands can be attributed to high mass and short-lived stars. 
As such any evolution in the SFR can be probed quite 
precisely, providing a suitably short redshift (i.e. temporal) 
baseline is chosen. Since large volumes are required to en- 
sure meaningful sampling statistics, large coverage areas are 
needed also — especially true for probing the local Universe. 
The local implied SFR is an extremely important normali- 
sation point for models that attempt to predict the cosmic 
star formation history (SFH), providing a value that is effec- 
tively independent of all but the most recent star formation. 
Since the local SFR can be measured with the greatest ac- 
curacy, it provides the main constraint to the normalisation 
of SFH. 

Much recent work has attempted to measure the local 
SFR, predominantly through calculating and integrat- 
ing out_Jhe_UV^ (IWvder et al.1 l2005l: iBudavari et al.1 
I2OO5I: 



IWvder et al 



20071) 



Ha 



or OH (|Gallego et al.1 I2OO2I I 



ISalim et al.l 

(jPerezGonzalez et al.l 
luminosity functions. However since these studies were pub- 
lished a huge amount of data has been released, increasing 
the number of galaxies with both measured UV fiuxes 
and redshifts. C onsidering purely F UV derived SFRs the 
mcrease is 15 X (IWvder et al.ll2005l \ although more recent 
work has used FUV data as part of a fitted SED derived 
SFR, makin g the data i ncrea se in this work slightly more 
than 1.3 X iSalim et aP l2007l V This work calculates local 
fiux densities in the redshift ranges 0.013 ^ z ^ 0.1 for the 
GALEX FUV and NUV bands, and calculates the local 
SFR based on the FUV fiux density. Throughout we use 



a cosmology of Ho =100kms~^, Qm — 0.3 and — 0.7. 
The work was carried out using a co mbination of publicly 
available GALEX MIS survey data (Morrissev et ajj l200'n ) 
and SDSS DR7 redshift data (Abazajian et al. 20ol). The 
GALEX MIS data are available from the MAST websit^B, 
and the SDSS DR7 is available from the CasJobs SQL 
websit^l- Since the FUV detecter on GALEX is no longer 
operational, the sample presented here is close to the largest 
that will be available in the near to mid future. 



2 GALEX-SDSS DATA PROCESSING 
2.1 Data Acquisition 

To generate our NUV and FUV luminosity function we 
require all of the potential object matches between the 
GALEX MIS and SDSS DR7. The first step was down- 
loading the full MIS catalogue for tiles with exposure times 
greater than 1,400 seconds and E{B — V) < 0.08, ensuring 
uniformly adequate depth and minimal Galactic dust cor- 
rections. To extract the overlap search region requi red for 
this work cone searches were made through Topcat (iTavlod 
2005), using the MIS tile centres and a search radius of 0.6°. 
Doing searches with Topcat proved to be a more efiicient 
technique than using the CasJobs service since the MIS tiles 
were spread out over such a complex survey geometry. We 
further cut this sample at an SDSS rpetro < 20.5mag, leaving 
a superset of 15M potential GALEX objects in the SDSS 
catalogue, and 6.4M primary photometry GALEX objects 
in the full MIS catalogue. Figure [1] shows the complex over- 
lapping coverage between the GALEX MIS and SDSS sur- 
veys. 

Due to the complex overlap between the GALEX MIS 
and SDSS DR7 surveys, calculating the area coverage is non 
trivial. To estimate the area we pixelated the common re- 
gion of a cylindrical projection and counted the number of 
joint GALEX-SDSS galaxies inside each pixel. Those pixels 
that had 8 immediate neighbours with 1 or more counts were 
considered to have complete coverage. The total number of 
objects contained within these 'complete' pixels was com- 
pared to the total number of objects over all pixels, and the 
total area of the 'complete' pixels was scaled by the same 
ratio. Whilst this process is sensitive to the pixel size chosen, 
there is a steady plateau of values for pixel sizes where the 
ratio between the number of 'complete' pixels and the num- 
ber of partial coverage pixels (1 or more counts) is highest. 
The pixel size where this statistic is highest should give the 
best estimate of the area, since the correction is smallest. 
The peak found for the GALEX-SDSS overlap data was for 
a pixel of size 15^ x 15^ at the equator, which gives 36.4% 
'complete' pixels. When the number count scaling is applied 
the total common GALEX-SDSS area is found to be 833.13 
deg\ 

To assess errors associated with our area estimate, sim- 
ulations were made of survey overlaps much like that used 
in this work — patches of contiguous overlapping GALEX 
size fields using the same number densities as found em- 
pirically. These were generated in areas ranging from single 
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Figure 1. Plot showing the overlap between GALEX MIS and SDSS DR7 redshifts. The total shared area is 833.13 sq. deg. 



GALEX field scales to 1000+ deg^, with a restriction that 
no contiguous block of GALEX fields can cover more than 
100 deg^- this is to ensure we investigate the worst case 
scenario of extremely patchy overlap. The results displayed 
marginal biasing towards overestimating the area, but errors 
are dominated by the noise in the estimation. The standard 
deviation in the area scatter is ^2.30% for areas below 50 
deg^, and drops to 0.34% for area between 800 and 900 deg^. 
This implies an area estimate of 833.13 ±2.83 deg^. In prac- 
tice the estimation should be better than this because the 
true GALEX-SDSS overlap contains much larger contiguous 
regions than allowed in the simulations, and the errors drop 
much faster for such structures than those implied by adding 
small areas and calculating the errors in quadrature. 



2.2 GALEX Flux Correction 

GALEX typically has FWHM imaging of 5^^ much greater 
than the 1.5^^ experienced with SDSS which we use for 
redshifts. To combat potential ambiguity in GALEX flux 
assignment a novel approach has been used based on the 
fact that the object centroid available for GALEX data are 
a flux moment based measurement. This means the centroid 
of the system is defined as the centre-of-mass (COM) in flux 
space. 

The COM in any system can be calculated by finding 
the COM of subsystems in order to simplify the problem, e.g. 
in the Earth-Moon system we can simplify the Earth and 
Moon as point sources and calculate the complete system 
COM accordingly. The same is true for calculating centroids 
in an image: in a symmetrical extended source (e.g. galaxy) 
the COM is at the peak of the light profile, and for a point 
source (e.g. star) it is also at the peak. If these two objects 
were incorrectly merged and the COM of the system was 
found to be half way between these two points then we would 
know that they each possess an equal amount of flux. For 
this work, if we only had access to the merged GALEX data 
then the problem would be intractable, but the associated 



SDSS data contains the centroid locations of all the likely 
component objects that were merged in GALEX. 

Since the object detection threshold used for GALEX 
will also be affected by significant de-blending issues when 
sources are close, we make a match for all object centroids 
(stars and galaxies) within the Petrosian rgo distance (the 
radius which encloses 90% of the r-band flux) for primary 
GALEX sources (this flag selects the best GALEX photom- 
etry in overlapping regions) and tag these as objects that 
potentially contribute some of the NUV / FUV flux from 
GALEX. Of the potential SDSS objects, 273k match am- 
biguously and 2.23M match unambiguously, so 83% of 
SDSS sources down to rpetro < 20.5mag have a GALEX 
object associated with it, of which 74% have an unambigu- 
ous 1-1 source match. At this stage the matching explicitly 
includes stars — this is important since often stars are the 
contributor to the UV flux and must be kept in the catalogue 
until after the flux has been redistributed. 

The flux for ambiguous targets is subsequently dis- 
tributed according to a few basic rules. Firstly, the SDSS 
object that matches closest to the GALEX centroid is des- 
ignated the primary match, and if no other matching object 
is within 2.5 magnitudes of the ^-band flux then all of the 
GALEX flux is assigned to the primary target (this proved 
to be a necessary addition to prevent large objects from be- 
ing shredded by spurious near-by SDSS detections). For this 
decision the ^-band was chosen due to its significant photo- 
metric fidelity compared to the i^-band data. Assuming some 
brighter sources remain as potential matches the next phase 
is to distribute flux according to the linear separation be- 
tween the GALEX and SDSS centroids, e.g. 

n 

i=l 



where is the separation between the GALEX source centre 
and the i^^ SDSS object, fr is the total GALEX flux, n is 
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Figure 2. Schematic showing the sort of ambiguous matches that 
occur between SDSS (red resolved galaxies) and GALEX (large 
blue unresolved detection). In practice such ambiguity is rare- the 
majority of matches are unambiguous, and there are only a few 
dozen GALEX sources with 4 or more SDSS matches. 



the number of SDSS objects that meet the matching criteria 
specified above and fi is the implied GALEX flux for the i^^ 
SDSS object. 

Figure [2] is a schematic of how the GALEX UV flux 
(the extended blue source) can be associated with multiple 
SDSS sources (the smaller red objects). The source clos- 
est to the GALEX centroid will take most of the flux- a 
reasonable assumption given that the SExtractor object 
coordinates are based on flux moments. The advantage of 
this technique compared to redistributing flux based on g- 
band flux, or possibly even SEDs, is that incorrect matches 
that made it through the initial g-band selection will barely 
affect the GALEX UV magnitude since they will be more 
likely to lie further from the centroid in relative terms. If 
prior magnitude knowledge is used to make flux redistribu- 
tion decisions this could introduce a significant bias into the 
derived colours. 

Using this simple technique of flux redistribution we 
find that the distribution of NUV — r colours for ambiguous 
sources moves very close to those calculated for unambigu- 
ous sources (see Figure [3]) . This is a very good indication 
that redistributing the flux in the manner described is both 
the right idea in theory, and has the desired empirical effect. 

A plausible concern of the approach taken is that large 
objects might become broken up spuriously (despite the g 
band selection designed to prevent this occurring). Consid- 
ering objects that make the final sample and have an ini- 
tial NUV ^ 17 (the brightest subset), - 40% of objects 
are associated with multiple SDSS objects that perhaps re- 
quire some amount of flux. Of galaxies corrected, the median 
NUV (FUV) magnitude shift is 0.31 (0.40), and in total 
this sample contains 98% (99.4%) of the total flux originally 
associated with it. For the full sample the corresponding fig- 
ures are NUV (FUV) magnitude shifts of 0.41 (0.47), and 
the whole sample contains 99.8% (99.8%) of the flux after 



flux redistribution. This would suggest that the flux for the 
brightest galaxies is not being reduced drastically compared 
to the full sample. 

A side effect of the flux redistribution process is the cre- 
ation of a very faint population of objects caused by erro- 
neous distant objects taking a small amount of a large flux. 
In these cases, the primary object flux will barely change 
but a notable number of objects will gain a negligible flux. 
Since these magnitudes are far beyond the apparent magni- 
tude limit used for producing the luminosity function, the 
effect is actually very minimal. 

Overall 79.2% of sources are unambiguous for our su- 
perset galaxy sample (SDSS objects flagged as galaxies), as 
deflned by the matching mechanism outlined above. When 
applying the flnal survey selection cuts (discussed in de- 
tail later) ^ 93% of the sample is unambiguously matched 
(varying slightly on the GALEX band). It should be high- 
lighted that applying this flux correction does not have a 
dramatic effect on the flnal results presented here, the in- 
tegrated fluxes remain almost the same. The main advan- 
tage of the approach advocated here would be for creating 
better quality SEDs for large samples of galaxies, and as 
such it has been applied t o the multi-wavelength GAMA 
( iDriver RobothanJ I2OIOI I survey which uses data from 
multiple PSF mis-matched sources. 

Whilst the ultimate result of the flux correction method 
described appears to be better than just doing a simple 
nearest-neighbour match, there are a number of obvious 
caveats. Firstly, it assumes that the proflle is circular. This 
will not always be the case, but for the majority of objects 
the 2D light proflle convolved with the GALEX PSF wifl be 
approximately circular. This does mean that flux could be 
incorrectly redistributed to an off-axis object that is within 
the rgo radius, when a simple match would not have done 
so. This is only likely to be a serious problem for the very 
largest objects rgo > 10", and these objects do not suffer 
unusually large amounts of flux loss in the redistribution 
process, suggesting the effect is not likely to be dominant. 
Another limitation is that this method does not account for 
colour in an explicit manner: galaxies which are very blue in 
g — r are likely to deserve more missing UV flux than galaxies 
which are very red in ^ — r. There are some indirect limits 
imposed by the requirement that the additional objects be 
within 2.5 mag of the primary match in the g-band, but oth- 
erwise there is no colour restriction in the implementation 
described. The use of colour (and by extension full SEDs) 
was considered, but there was a desire to keep the algorithm 
as simple as possible. Figure [3] suggests that the red popula- 
tion is the one most overtly affected by flux loss in GALEX, 
and in fact more flux should be shifted from the blue to red 
populations. Despite the limitations discussed, flux shifting 
appears to produce a better fldelity of data than simply 
matching via nearest-neighbour. 

2.3 Completeness 

To calculate the completeness function required to gen- 
erate galaxy weights the SDSS star galaxy classiflcation 
was adopted to identify potential targets that could have 
been selected by virtue of their r-band luminosity (rpetro < 
17.77mag). GALEX target completeness was identifled as 
a function of rpetro in the early and late type populations 
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Figure 3. Effect of the flux correction applied to ambiguous matching sources between GALEX and SDSS. The plots show how the 
correction affects the FUV — r colour (left) and NUV — r colour (right) for galaxies with redshifts in our sample. Black line shows 
the distribution for unambiguously matched SDSS-GALEX sources, the dashed red line shows the ambiguous sources before correction, 
and the solid red line after correction. It is particularly noticeable in the NUV that the post corrected distribution is much closer to 
the unambiguous source matched distribution. The bimodality also strengthens significantly when the correction is applied, particularly 
evident for NUV — r colours. 
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Figure 4. NUV — r as a function of redshift. The blue line rep- 
resents the peak in the bimodality for blue (late- type) galaxies, 
the red line represents the peak in the bimodality for red (early- 
type) galaxies and the green line represents the minima between 
the two peaks. The linear trend seen for blue (late- type) galaxies 
over this redshift regime is just a selection effect caused by sam- 
pling lower luminosity (more numerous) objects at lower redshifts 
for a fixed blue colour. 



separately. To split the populations the NUV — r colour 
bimodality was studied as a function of redshift, and the 
minimum in the bimodality was used to select late-type star 



forming populations and redder colours to select eary-type 
quiescent populations (see Figure 2]). The dichotomy is sig- 
nificant in terms of completeness, with late- type galaxies 
possessing a higher completeness by virtue of their easily 
identifiable emission features aiding redshift measurements. 
Over the vast majority of galaxies the completeness is quite 
flat as both a function of NUV or r, where the average 
completeness for early- type galaxies is ^ 78%, for late- type 
galaxies it is ^ 84% and globally it is ^ 82%. The inverse of 
the completeness as a function of rpetro and galaxy type is 
used to generate weights for each galaxy. Assuming a con- 
stant a magnitude dependent weight does not have a large 
impact on any of the results presented here. 

2.4 Photometric Corrections 

The GALEX photometry acquired from MIS required Galac- 
tic extinction corrections and K-corrections. For the for- 
mer E{B — V) values calculated using the Schlegel (1998) 
dust maps. The reddening factor for GALEX bands is 
Afuv/E{B -V) = 8.24 mag and Anuv/E{B - V) = 8.2 
mag {NUV has a sligh t spectral dependen ce, but we use the 
constant values used in 'Wvder et al.l l2Q07l ). which due to our 
low extinction plate requirement means the mean correction 
is less than 0.34 mag in both bands. 

The k-corrections appl i ed wer e created using kcorrect- 
V4.2 (jBlanton Roweid l2Q0/t ). which fits the dust 
corrected photometric SED for each galaxy (FUV^ron, 

NU Vhron, UpetTO, 5'petro, Tpetro, Rpetro, 2;petro WCrC the bauds 

used in this work, and should be presumed throughout) to a 
library of templates and calculates the change in flux mea- 
surement due to the the filters shifting and stretching com- 
pared to observing the same galaxy at z = 0. Figure [5] shows 
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Figure 5. Comparison of FUV — r and NUV — r colours before (upper panels) and after (lower panels) k-corrections. The 1,000 points 
in the least dense part of the parameter space are plotted, for the rest the intensity of colour is a guide to the local density of points. 
The significant observation is how the bimodality strengthens and the redshift dependence disappears post correction. 



the effect of applying the k-correction to the FUV and NUV 
bands as a function of redshift. The most noticeable effect 
is that the early- type/ late- type bimodality most evident in 
NUV — r becomes more distinct. The NUV — r colour for 
redder (early- type) galaxies is much more distinct, and ap- 
pears to remain very flat as a function of redshift. The bluer 
(late-type) galaxies in the NUV — r plot appear to show 
strong evidence for evolution in colour over the baseline of 
3.4 Gyr (out to z = 0.3) shown in this plot, but this effect 
is actually a selection effect caused by the fixed r and NUV 
limits. Crudely speaking an object with a NUV — r ^ 2 with 
r = 17.77 (i.e. at the SDSS limit) will have Mnuv=-13.21 
at z — 0.013, whereas at z — 0.3 the same magnitudes cor- 
respond to Mnuv=-20.41. Assuming a typical distribution 
of absolute magnitudes, we would expect fewer galaxies to 
be present at z — 0.3 for the same NUV — r colour. We do 
not attempt to correct for evolution in the sample, the cor- 



rections being small over the ^1 Gyr baseline investigated 
here. 

2.5 Sample Limits 

Choosing a sensible depth for the survey, and mapping 
this onto an observable absolute magnitude limit, is impor- 
tant. The process is made more complicated by virtue of 
the redshift catalogue survey (SDSS) is limited by rpetro < 
17.77 mag, whilst GALEX is photometrically limited in the 
FUV and NUV bands. The colour magnitude diagrams for 
FUV — r versus FUV and NUV — r versus FUV are shown 
in Figure [6l The assumed limits are 23.0 mag for both FUV 
and NUV, this i s brighter t han t he 23.8 mag and 23.6 mag 
depths given in IXu et al.1 (|2005l ) for FUV and NUV so 
should be a conservative estimate. 

Due to the limited depth of the SDSS redshift survey 
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Figure 6. Comparison of FUV and NUV survey depths. The effective rpetro hmit is 17.77 mag throughout, and the assumed maximum 
GALEX hmits are 23.0 for FUV and 23.0 for NUV (these are the dotted hues drawn). The floating survey hmit in each GALEX band 
is that imphed by the sohd black hue, where the GALEX colour (FUV — r or NUV — r) and the GALEX apparent magnitude deflne 
the survey limit. 



(median redshift is z ^ 0.1), and the desire to sample a 
minimal amount of evolution so as to obtain flux densities 
for the most local Universe, we used a redshift sample span- 
ning 0.013 ^ z ^ 0.1, where the low redshift cut is used to 
avoid nearby galaxies that require signiflcant local flow cor- 
rections, and the hi gher redshift limit is the same as used in 
IWvder et al.l (|2005l l . This makes comparison more straight- 
forward, and keeps the measured UV flux within a 1 Gyr 
baseline. 



3 RESULTS 

The Luminosity Distributions (LDs) for GALEX were gen- 
er ated using a standa rd l/Vmax method (in a manner similar 
to lWvder et al.l [20071 ) for 10 subsets of the sample split by 
NUV — r colour. The subsets were determined by calcu- 
lating the even quantiles of the data. To minimise colour 
bias within each subset of the data we calculated the 10*^ 
percentile of each subset and use this as our effective magni- 
tude limit for NUV /FUV , leaving the r-band limit fixed at 
Rpetro ^ 17.77 mag. Attempts were made to split the sample 
up into either simply late/early type and late/central/early 
type, but the colour bias was such that the resulting LPs 
were heavily distorted. Using the colour binning specified 
was a good compromise between data size for each subset 
and quality of the resulting LD (i.e. by visual inspection it 
appeared physical). Such an approach also means we keep 
the maximum possible number of galaxies when calculating 
the total integrated luminosity density of the local Universe. 

The consequence of dividing the sample like this is that 
we are not left with a simple blue (late-type) and red (early- 
type) LDs. To determine appropriate luminosity densities 
for the full sample, and the early/late-type populations, 
stacking of these distributions must be done. Since we are 
interested in the integrated light the best approach is to fit a 



parametric function to the data subsets and to integrate un- 
der this function beyond the depth we can observe to. Since 
we already know the fraction of early/late-type galaxies be- 
longing to each colour subset, we can share the integrated 
light between the two populations to calculate the luminos- 
ity density for the early-late-type and composite populations 
separately. 



3.1 Luminosity Functions 

To these empirical LDs we fitted the analytic Schechter 
Luminosity Function (LF), a 3 parameter function that is 
widely used in the literature, and allows for s imple compar- 
isons between models to be made (|SchechterllT976 ): 



(3) 



A more complicated form of this function can be defined 
in terms of magnitudes (M* and M) in place of luminosities 
(L* and L). The 3 parameters that we fit are M*, which 
for most distributions represents the magnitude of the ob- 
ject that provides the most overall fiux; 0*, approximately 
the comoving space density of galaxies with the magnitude 
M* and acts as a normalisation constant; and a, a power 
law term that describes the faint end slope of the LF. The 
combination of all 3 parameters allow us to calculate the en- 
ergy output per unit volume. Even though the faint end will 
truncate eventually, the fiux output is usually dominated 
by galaxies near M* , so the important consideration is that 
we sample beyond M* to ensure this is well measured. In 
practice 1 mag is the minimum beyond M* that should be 
sampled. To fit the Schechter function we find the minimum 
fit for each distribution, where the errors are calculated 
using the covariance matrix produced during the fit. 

The empirical LDs created via l/Vmax and the best fit 
Schechter LFs can be seen in Figure[71for the 0.013 ^ z ^ 0.1 
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Figure 7. Comparison of FUV and NUV luminosity functions for early-type, late-type and composite populations of galaxies. The 
blue-green-red colour scaling maps maps onto the bluest (lowest) FUV — r (NUV — r) to the reddest (largest) FUV — r (NUV — r). 



sample. The full table of Schechter parameters and derived 
values is presented in Tables [T] and (2] The median x^/^ for 
each fit is 1.07/1.02 for NUV and FUV data respectively, 
demonstrating that the subset divisions used are appropri- 
ate since the data are neither under or over- fitted. Figure [9] 
shows the 1, 2 and 3a error contours for the colour subsets 
used in this work. A dramatic feature of the LFs is that they 
maintain a relatively constant M* for the first five colour 

bins for both FUV and NUV (M^uv 17.2 mag and 

^Nuv —17.7 mag) but with a rapidly increasing (-1.5 - 
-0.7) a. The value of a reaches a maximum value for the 
jth i^th (3qJq^j. subsets combined with a rapidly faintening 
M*, before a appreciably steepens again for the final few 
reddest bins. 

The transitions in M* and a are statistically 
significant — the major parameter shifts actually occur or- 
thogonally to the main direction of M* and a degeneracy 
(the vector of ellipse elongation seen in Figure (9]) . Together 
with the shifts in 0*, these fits mean the integrated flux 
varies hugely as a function of FUV — r of NUV — r colour. 
The immediate observation is that the bluest colour bins 
contain the vast majority of the integrated flux, this is par- 
ticularly clear in Figure [8] which shows the fractional contri- 
bution of each colour subset to the total integrated light. The 
early/late-type contributions are also shown, where 97.9% 
(95.5%) of FUV (NUV) light is found in late-type galaxies 
in the local Universe. 

The FUV 3nd NUV LFs pres ented here have similari- 
ties to those in lWvder et al.l (|2007[ ). Figure 19 in that paper 
can be compared almost directly with the NUV LFs in Fig- 
ure [71 Both demonstrate a significant faintening of M* and 
a flattening of a as the subsets transition from blue to red 
populations. IWvder et al] (l2007l l use regular colour subsets 
to divide the populations, whereas here the number of galax- 
ies used to construct each LF is designed to be as similar as 
possible, as such the most extreme red and blue populations 
investigated in that work are not present here. 




NUV-r/FUV-r 



Figure 8. Fraction of total luminosity density contained for each 
colour subset. Dashed and dotted lines show the late/early- type 
contribution separately, where, as should be expected, the late- 
type population dominates the bluer colours. Overall, the vast 
fraction of integrated UV light (FUV and NUV) comes from the 
bluer subsets/ late- type populations. 



3.2 Composite Luminosity Function 

To allow for easy comparisons to further work the full com- 
posite galaxy densities and errors have been calculated using 
0.25 mag bins spanning -20.25 to -14.00 for both FUV and 
NUV data, using the 10 colour split populations shown in 
Figure [71 These are presented in the Appendix, and include 
summed versions of the early-type and late-type population 
samples for each magnitude bin. The preference is that these 
empirical values should be used in any comparison, but for 
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Figure 9. Comparison of FUV (left) and NUV (right) luminosity functions for the 10 different colour subsets given in Tables [Hand El 
The blue-green-red colour scaling maps maps onto the bluest (lowest) FUV — r (NUV — r) to the reddest (largest) FUV — r (NUV — r). 



completeness the best composite FUV Schechter function 
has 0* = 0.0165 ± 0.0007/i^Mpc-^mag-\ M* = -17.031 ± 
0.030 and a = - 1. 167 =b 0.025 and the best composite NUV 



M* 



-17.641 ±0.025 and a 



-1.200 ±0.019. 



3.3 Luminosity Densities 

The most interesting derived components, given a well be- 
haved LF that can be parameterised accurately with a 
Schechter fit to the distribution, is the total integrated 
light output of galaxies per unit volume and the implied 
star formation density (SFR) of the Univers e. Following the 
approach taken in many s t udies (e .g. M adau et al. 1 19961: 
Baldrv Glazebrookl l2003l : IPriver et all l2008l : iHill et al.1 
20101 ). we can integrate the implied Schechter function out 
to infinity using the F function: 



fu = h(10 



-(M*-2.51ogio(0*r(2+a)))/2.5 



)4.345 X 10^ 



(4) 



This will give the Universal flux density in units of ergs 
s~^ Hz~^ Mpc~^, and because of the cancelation between 
absolute magnitude (which depends on distance squared) 
and density (distance cubed) is only dependent on a sin- 
gle factor of /i, assuming otherwise standard cosmology. As 
discussed above, the influence of the faint-end is relatively 
minor for typical a values. If the empirical LF extends 4 
magnitudes faintwards of M* then integrating out to in- 
finitely faint fluxes rather than using the limit implied by 
the data only increases the flux density by 2.5% for a = —1 
and 13.8% for a = —1.3. At worst, this effect is of similar 
importance as the cosmic variance in the sample. The role 
of cosmic variance (CV) for different surveys, and in partic- 
ular different survey geometries, is discussed in some detail 
in Driver & Robotham et al. (2010). Taking our redshift 
range of 0.013 ^ 2; ^ 0.1 and the area coverage of 833.13° 
we have a total survey volume of 2.112 Gpc^h"^. Assuming 
the worst case scenario of our survey having equal dimen- 
sions the formula presented in Driver & Robotham (2010) 



predicts a CV^ 12.6%, and no lower than 12.0% assuming 
geometry similar to that actual used. All densities, and val- 
ues derived from densities, can be considered to possess this 
extra form of error. 

The derived parameters, fi^ and the SFR (deter- 
mined using (Kennicutt 1998) SFRFuv{MQyT~^)= 1.4 x 
10"^^L^(ergs s~^ Hz~M ), are largely consistent with those 



found by (Wvder et al.1 l2005l ) and (jBudavari et al.ll2005l l. 

the appropriate conversions to the units used in this paper 
are in Tables [1] and (2] The flux density in the FUV agrees 
with in implied errors for both the low r edshift (Wy de r et al.l 
I2OO5I ) sample and the higher redshift (jBudavari et al.ll2005l l . 
There is more of a discrepancy for the NUV where we find 
f^uy = 9.706 ± 0.541 X 10^^ (h ergs s"^ Hz"^ Mpc-^) 
and Wvd er et al.1 (|2005l ) find /nuv = 7.52 ± 2.16 x 10^^ 
(h ergs s~^ Hz~^ Mpc~^). This is only slightly larger than 
the formal errors calculated from integratin g the LFs. A 
possible cause for this difference is that the IWvder et al. 
(2005) redshift sample was constructed from matching the 
overlap between the GALEX MIS and the Two Deg ree 
Galaxy Redshift Survey (2dFGRS; IColless et all l200li i. a 
co mpletely separa t e volu me to the one considered here and 
in lBudavari et aL I (l2005h . This means differences arising due 
to CV effects are expected to be more significant since there 
is no common volume probed in the two studies, an effect 
that is discussed in detail la ter. The NUV flux density esti- 
mate in I Wvder et al.l (|2007f ) is much closer to that presented 
here: fNuv = 8.83 ± 0.61 x 10^^ (h ergs s"^ Hz"^ Mpc"^). 
This later paper uses SDSS data also, so a convergence in 
results should be expected since a large fraction of the data 
will be the same between the two analyses. 



3.4 Dust Corrections 

All the flux density numbers presented so far are uncor- 
rected for dust attenuation within the host galaxy (Galactic 
extinction has already been accounted for). Much work has 
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FUV - r 


0* 


M* — 5 logiQ h 




fi, X 10^^ 








(mag) 




(h ergs s~^ Hz~^ Mpc~"^) 




2.215 


0.00345 ± 0.00046 


-17.082 ± 0.086 


-1.436 ± 0.058 


1.607 ± 0.288 


1.835 


2.607 


0.00266 ± 0.00026 


-17.111 ± 0.072 


-1.193 ± 0.052 


0.936 ± 0.114 


1.604 


2.881 


0.00236 ± 0.00018 


-17.064 ± 0.066 


-1.005 ± 0.052 


0.687 ± 0.073 


0.506 


3.134 


0.00187 ± 0.00014 


-17.146 ± 0.068 


-0.938 ± 0.055 


0.569 ± 0.059 


0.963 


3.390 


0.00207 ± 0.00011 


-16.862 ± 0.051 


-0.695 ± 0.060 


0.448 ± 0.031 


2.424 


3.683 


0.00219 ± 0.00009 


-16.526 ± 0.054 


-0.446 ± 0.070 


0.345 ± 0.023 


2.207 


4.040 


0.00217 ± 0.00008 


-16.253 ± 0.052 


-0.402 ± 0.061 


0.268 ± 0.017 


1.080 


4.532 


0.00210 ± 0.00009 


-15.843 ± 0.061 


-0.405 ± 0.069 


0.177 ± 0.013 


0.693 


5.207 


0.00224 ± 0.00011 


-15.205 ± 0.061 


-0.534 ± 0.065 


0.104 ± 0.008 


0.637 


6.562 


0.00339 ± 0.00030 


-14.497 ± 0.074 


-1.021 ± 0.069 


0.094 ± 0.011 


0.934 


Total 








5.234 ± 0.327 





Table 1. Table of Schechter Fitting Parameters for current work for FUV data split by FUV — r colour. The final row is the flux density 
for the composite population of all the above FUV — r subsets. 



NUV - r 


0* 


M* — 5 log]^Q h 




fi, X 10^5 






(h-^Mpc^-^mag"^) 


(mag) 




{h ergs s~^ Hz~^ Mpc~"^) 




1.868 


0.00279 ± 0.00039 


-17.719 ± 0.088 


-1.495 ± 0.051 


2.608 ± 0.472 


2.046 


2.203 


0.00291 ± 0.00026 


-17.613 ± 0.068 


-1.158 ± 0.047 


1.574 ± 0.181 


0.973 


2.444 


0.00206 ± 0.00017 


-17.773 ± 0.068 


-1.120 ± 0.045 


1.248 ± 0.132 


0.963 


2.662 


0.00189 ± 0.00014 


-17.767 ± 0.066 


-0.976 ± 0.049 


1.036 ± 0.105 


0.970 


2.881 


0.00222 ± 0.00011 


-17.500 ± 0.051 


-0.688 ± 0.057 


0.862 ± 0.062 


2.050 


3.134 


0.00225 ± 0.00010 


-17.290 ± 0.051 


-0.571 ± 0.055 


0.715 ± 0.046 


1.965 


3.453 


0.00234 ± 0.00010 


-16.996 ± 0.056 


-0.477 ± 0.060 


0.568 ± 0.038 


1.181 


3.882 


0.00252 ± 0.00009 


-16.531 ± 0.048 


-0.377 ± 0.059 


0.402 ± 0.023 


1.742 


4.475 


0.00258 ± 0.00013 


-16.235 ± 0.054 


-0.723 ± 0.047 


0.315 ± 0.024 


0.887 


5.548 


0.00828 ± 0.00023 


-15.189 ± 0.034 


-0.500 ± 0.034 


0.379 ± 0.015 


0.808 


Total 








9.706 ± 0.541 





Table 2. Table of Schechter Fitting Parameters for current work for NUV data split by NUV — r colour. The flnal row is the flux 
density for the composite population of all the above NUV — r subsets. 



been invested into the attenuation caused by dust, and the 
effect has been investigjated in detail in rec ent papers by 
iTuffs et alJ (|2004l l: IPriver et alJ (|2007l . l2008l ). The implied 
typical (so not taking into account specific galaxy geometry, 
but the mean effect of dust for an integrated population) cor- 
rection for the FUV is reported to be a f actor of 4.32, and fo r 
the NUV a factor of 4.11, as taken from [Driver et al.1 (|2008[ ). 
This would suggest the unobscured flux density for our sam- 
ple is fr^-Fuv = 22.61 ±1.41 xlO^^/iergss-^ Hz"^ Mpc-^ 
f^_Nuv = 39.89 ± 2.22 xlO^^h ergs s"^ Hz"^ Mpc-^ 

A bette r estimate can be made by only applying the 
iDriver et all ([2008. ) corrections to the late- type selection, 
since we would typically expect little dust in our early-type 
sample. This gives a dust corrected 0.013 ^ z ^ 0.1 late- 
type flux density of fu-Fuv = 22.13 ± 1.41 xlO^^/i ergs s"^ 
Hz-^ Mpc-^ fu-Nuv = 38.10±2.22 xlO^^h ergs s"^ Hz"^ 
Mpc~^. If these dust corrected values are added to the non- 
dust corrected early-type sample we can produce a final flux 
density for the Universe, with proper consideration made 
for the average dust attenuation. It is worth highlighting 
that the late-type galaxy sample dominates the integrated 
flux, thus the effect of only applying dust correction to the 
early-type galaxies changes the flux densities by only a few 
percent. Table [3] presents the final integrated flux densities 
for the local Universe for FUV and NUV, and the implied 
SFR. 



3.5 Low-z SFR Comparison 

The SFRs quoted in Table [3] are somewhat higher than 
typically reported for UV inferred values taken straight 
from the literature, but they are consistent with the 
SFR of 0.029 + 0.008 - 0.005 h Moyr-^Mpc"^ given by 



'PerezGonzalez et al. 1 (12003"). It is clear from this table that 
with data currently available we are dominated hugely in any 
estimate of the local SFR by cosmic variance. Since the light 
output in the GALEX UV bands is dominated almost solely 
by short lived stars (age less than 100 Myr) it is not sensible 
to probe over deeper baselines in z where age resolution will 
be washed out. For measuring the cosmic star formation his- 
tory (CSFH) this is a particular concer n since we expect it 
to drop so steeply from z = 1 to z = ([Hopkins &: BeacomI 
2006). The limit on measuring the SFR at the 'present' day 
will always be the cosmic variance of the local Universe, a 
rough prediction from the Driver & Robotham CV calcula- 
tor suggests that within z = 0.08 (the redshift limit out to 1 
Gyr using Ho = 71kms"VMpc, = 0.27 and Qa = 0.73) 
the local SFR cannot be constrained to better than 4.3% 
even if we could measure the UV flux density in the whole 
sky, a survey area 50 times larger than that used for this 
work. 

In an eflFort to consolidate the various low redshift mea- 
surements of the SFR we conducted an up-to-date literature 
search, normalising the quoted SFRs to our Hq = lOOkms"^ , 
Qm = 0.3 and Qa = 0.7 cosmology. Further to this, where 
UV SFRs were quoted, the Driver et al. (2008) dust correc- 
tions were applied to the data in order to treat the data 
as consistently as possible. Where possible, the volume of 
the various surveys were determined based on quoted area 
coverage and redshift range, allowing us to calculate the im- 
plied CV of each value using the same Driver & Robotham 
(2010) equation as discussed above. If sharp redshift bound- 
aries were not used, then the median redshift of the survey 
was used to determine an approximate volume. These addi- 
tional CV error estimates were added in quadrature to the 
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z-range 
z-range 


/i. X 10^5 (original) 
h ergs s~^ Hz~^ Mpc""^ 


fi, X 10^^ (corrected, Driver 2008) 
h ergs s"-*^ Hz""*^ Mpc~^ 


SFR 

h MQyr~-'-Mpc~'^ 


FUV composite 
FUV late 
FUV early 


0.013-0.1 
0.013-0.1 
0.013-0.1 


5.234 ± 0.327 ± 0.660 
5.123 ± 0.327 ± 0.646 
0.111 ± 0.010 ± 0.014 


22.24 ± 1.41 ± 2.80 
22.13 ± 1.41 ± 2.79 


0.0312 ± 0.0020 ± 0.0040 
0.0310 ± 0.0020 ± 0.0039 
0.00016 ± 0.00001 ± 0.00002 


NUV composite 
NUV late 
NUV early 


0.013-0.1 
0.013-0.1 
0.013-0.1 


9.706 ± 0.541 ± 1.170 
9.270 ± 0.540 ± 1.169 
0.436 ± 0.016 ± 0.055 


38.54 ± 2.22 ± 4.81 
38.10 ± 2.22 ± 4.80 




FUV Wyder (2005) 
NUV Wyder (2005) 


0-0.1 
0-0.1 


5.17 ± 0.98 ± 1.78 
7.52 ± 2.16 ± 2.59 


22.33 ± 4.23 ± 7.70 
30.91 ± 6.21 ± 10.65 


0.0313 ± 0.0059 ± 0.0108 


FUV Budavari (2005) 
NUV Budavari (2005) 


0.07-0.13 
0.07-0.13 


4.54 ± 1.17 ± 1.38 
7.25 ± 1.72 ± 2.20 


19.61 ± 5.05 ± 5.96 
29.80 ± 7.07 ± 9.04 


0.0275 ± 0.0071 ± 0.0083 


NUV Wyder (2007) 


0.01-0.25 


8.83 ± 0.61 ± 0.32 


36.29 ± 2.51 ± 1.31 





Table 3. Table of dust corrected flux densities, and the derived iKennicuttl (|l998l ) SFRs. The flrst part of the errors contain the covariance 
errors carried through from various sources, and the second part include a 12.6% error expected due to Cosmic Variance (see Driver 
& Robotham 2010) and a 0.36% error in our area estimation (for this work), a 34.5% CV error for Wyder 2005, 30.4% CV error for 
Budavari 2005 and 3.6% CV error for Wyder 2007. The CV error dominates over all other sources for all but Wyder 2007. 



random errors stated for each SFR, allowing us to compare 
how unexpected any deviations in SFRs truly are. 

Table 2] presents the relevant figures for these differ- 
ent SFR measurements, along with our new CV values. 
The table also explains what assumptions were made to 
obtain the SFR, i.e. what type of dust correction was ap- 
plied (Driver, Calzetti or non-specific Balmer decrement cor- 
rection) and what co nversion to SFR was used: jKennicuttl 
([l998 ) or PEGASE (iFioc Rocca-Volmerangdl 19971 ). Fig^ 
ure [To] compares the SFRs as a function of redshift, sep- 
arating out the quoted random errors found in the liter- 
ature (small error bars) to the larger errors implied by 
adding the worst case scenario CV in quadrature (large error 
bars). The sour c es for t his compendiu r a of l o w-z SFRs are 
iBudavari et al. (I2OO5II. iGallego et al.1 (I2002T) Hanish et al. 
(2006'), "PerezGonzalez et al." (2003), "Salim et al." (2007), 
,Schriminovich et al. (2005), Sullivan et al. (2000) and 
[Wvder et al.1 (|2005h . These derive SFRs from utilising data 
from the following survey s: 2dFGRS, FOCA (|Mimard et al.1 
I1992I). GALEX . HIPASS (Mever et al.ll2004l l. SDSS, SINGG 
(Hanis h et "aD [20061 ) an d the Univ ersidad Complutense de 
Madrid Survey f UCMS [Zamorano et al. 1994). 

Figure [10] suggests that once CV is taken into account 
most of the values for the local SFR found in the literature 
agree well with the new, much more tightly defined, SFR 
presented in this work. It should be highlighted that the CV 
quoted is the worst case scenario CV, calculated by convert- 
ing the quoted volumes to cubes. Any geometry different to 
this will lower the CV, and actual values will often be a few 
percent smaller than stated here (see Driver & Robotham 
2010 for a detailed discussion of these effects). Another sub- 
tle effect is that some of the volumes used in various studies 
overlap to some degree, lowering their CV with respect to 
each other. Regardless, the CV is in practice always the 
dominant source of discrepancy, and the spread of values 
found in the literature can be explained by accounting for 
this source of error. 

The lowest redshift SINGG derived Ha and UCM de- 
rived OH SFRs do appear to be significantly low. Since these 
values probe significantly different redshift ranges (as indi- 
cated by the small error bars in Figure [TO]) , this could be in- 
terpreted as a real drop in cosmic SFR over the ^ SOOMyrs 
that covers the median redshift of our low redshift sample 
and median redshift of the two studies highlighted. 



Perhaps the only significant difference between SFR in- 
dicators comes from using OH as opposed to UV or Ha. 
The OH SFR values are generally lower than all three in- 
dicators, and using the same surveys (UCMS and FOCA) 
they are lower than values recovered via Ha. For the UCMS 
derived SFRs the difference is well outside the random er- 
rors quoted (the appropriate errors to compare since the two 
measures share the same volume and suffer CV in the same 
direction) . 



4 CONCLUSIONS 

In this paper we have a calculated a number of flux densities 
for both dust-corrected and non-dust-corrected by integrat- 
ing out to infinity various Schechter luminosity functions. 
We find the data are much better represented by dividing 
the populations into intense star forming populations and 
less active galaxies via their NUV — r colour. This possesses 
a very strong bimodality, and allows dust corrections to be 
applied more appropriately to just the population of galax- 
ies that are rapidly producing new stars and by inference 
will have the most dust attenuation. 

Whilst flux densities, and the SFR they imply, have 
been calculated for a variety of stages of analysis, the fi- 
nal local SFR found for the low redshift 0.013 ^ z ^ 0.1 
Universe is 0.0312 ± 0.0045 h Moyr^^Mpc"^. The majority 
of the quoted uncertainty is due to cosmic variance (often 
called sample variance) , which will dominate any calculation 
of the local SFR since it is strongly volume limited. 

Various other low redshift SFRs were obtained from 
the literature and compared to values found in this work. 
Once CV is properly accounted for the values generally agree 
within errors, with the possible exception of the lowest red- 
shift values — these appear to be significantly on the low 
side. Also, the OH SFR indicator appears to predict gener- 
ally lower values for the local SFR. 
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Figure 10. Comparison of low redshift SFR indicators taken from the literature when factoring purely random errors (left side plot) or 
considering CV also (right side plot). Colour shows indicator type (UV, Oil or Ha), and point style shows indicator source. Black boxes 
represent the results of this work (small error bars are the solid line, large error bars are the dotted line) On the x-axis, small error bars 
show the redshift at which 15% and 85% of the survey volume is reached, larger error bars show the extremities of the survey z-range. 
On the y-axis, error bars show the random errors quoted in the literature for the left side plot, error bars show the error implied when 
the worst case scenario CV error is added in quadrature for the right side plot. Exact values can be found in Tabled 



z-range SFR Area Volume X 10 CV Indicator Dust SFR Source Survey Source 

h MQyr^^Mpc"-^ (h^Mpc~^) % 



This work 


0.013-0.1 


0.0312 ± 0.0020 


833.13 


2,121 (exact) 


12.6 


UV 


Driver 


Kennicutt 


GALEX-SDSS 


Budavari 2005 


0.07-0.13 


0.0276 ± 0.0069 


43.9 


202 (exact) 


30.4 


UV 


Driver 


Kennicutt 


GALEX-SDSS 


Gallego 2002 


0-0.045 


0.014 - 0.005 + 0.005 


471.4 


114 (exact) 


37.6 


on 


Implied Balmer 


Kennicutt 


UCMS-HIPASS 


Hanish 2006 


0-0.042 


0.016 - 0.0041 + 0.0028 


23,504 


4,629 (scaled) 


9.4 


Halpha 


Helmboldt 


Kennicutt 


SINGG 


Perez-Gonzalez 2003 


0-0.045 


0.041 + 0.008 - 0.007 


195.0 


47.1 (scaled) 


56.2 


Halpha 


Implied Balmer 


Kennicutt 


UCMS 


Salim 2007 


0.005-0.22 


0.0261 - 0.0006 + 0.0021 


645 


16,046 (exact) 


6.1 


SED 


Chariot & Fall 


Bruzual & Chariot 


GALEX-SDSS 


Schiminovich 2005 


0-0.1 


0.0299 - 0.0048 + 0.0009 


3.141 


8.02 (exact) 


95.2 


UV 


Driver 


Kennicutt 


GALEX- VVDS 


Sullivan 2000 


0.01-0.3 


0.0573 - 0.0058 + 0.0070 


2.2 


131 (exact) 


35.8 


UV 


Calzetti 


Pegase 


FOCA 


Sullivan 2000 


0.01-0.3 


0.0256 - 0.0033 + 0.0038 


2.2 


131 (exact) 


35.8 


Halpha 


Calzetti 


Kennicutt 


FOCA 


Sullivan 2000 


0.01-0.3 


0.0223 - 0.0024 + 0.0027 


2.2 


131 (exact) 
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145 (exact) 


34.5 
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Table 4. Table comparing different SFR measures at 2; ^ 0.3 taken from the literature. Survey source shows the main combination of 
surveys required to generate the SFR. For most studies the secondary source supplies r edshifts to allow LF measur ements, in the case of 
Gallego (2002) the seconda ry source (HIPASS ) is used to calibrate the normalisation. Schriminovich et al.l (|2005l ) is effectively derived 
from the measurements of (jWvder et al.llioosh . so should not be considered a fully independent measure of SFR. 
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Table 5. Summed galaxy densities for FUV. Absolute magnitudes are for the bin centres. The columns show the summed densities for 
early, late and all galaxies. 
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Table 6. Summed galaxy densities for NUV. Absolute magnitudes are for the bin centres. The columns show the summed densities for 
early, late and all galaxies. 
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